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Abstract 

We discuss the causal diagrams of static and spherically symmetric bigravity vacuum 
solutions, with interacting metrics / and g. Such solutions can be classified into type I (or 
"non-diagonal") and type II (or "diagonal"). The general solution of type I is known, and 
leads to metrics / and g in the Schwarzschild-(Anti)de Sitter family. The two metrics are 
not always diagonalizable in the same coordinate system, and the light-cone structure of 
both metrics can be quite different. In spite of this, we find that causality is preserved, 
in the sense that closed time-like curves cannot be pieced together from geodesies of both 
metrics. We propose maximal extensions of Type I bigravity solutions, where geodesies 
of both metrics do not stop unless a curvature singularity is encountered. Such maximal 
extensions can contain several copies (or even an infinite number of them) of the maximally 
extended " individual" geometries associated to / and g separately. Generically, we find that 
the maximal extensions of bigravity solutions are not globally hyperbolic, even in cases when 
the individual geometries are. The general solution of type II has not been given in closed 
form. We discuss a subclass where g is an arbitrary solution of Einstein's equations with a 
cosmological constant, and we find that in this case the only solutions are such that f (x g 
(with trivial causal structure). 
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1 Introduction 



Recently, massive gravity, and other infrared modifications of General Relativity, have been the 
subject of numerous investigationsjTJ|21ISllllE], largely motivated by the enigma of dark energy 
in IH]. A linearized ghost-free theory of massive gravity can be constructed by adding a suitable 
quadratic mass term to the linearized Einstein Lagrangian. This is the so-called Pauli-Fierz 
theory (PF) P|. Non-linear and generally covariant extensions can easily be formulated, at the 
cost of introducing new fields. The simplest possibility is to consider the theory of two metrics, / 
and (7, with standard kinetic terms, where the mass terms arise from non-derivative interactions 
between / and g. This is known as bigravity, f — g theory, or strong gravity, and was first 
considered in the seventies, in the context of the strong interactions Massive gravity poses 

certain difficulties. First of all, the linearized theory does not agree with observations. The 
PF Lagrangian propagates a scalar degree of freedom (in addition to two tensor and two vector 
polarizations) which couples to matter with gravitational strength. The coupling persists in the 
limit as the graviton mass vanishes, causing the so-called vDVZ discontinuity [T5|. It has been 
suggested that the discontinuity may disappear due to non-linear effects in the sense that 
solutions of massless gravity can be recovered non-perturbatively. This mechanism may indeed 
work in more complicated infrared modifications of gravity, such as the theories discussed in 
although its implementation in the context of bigravity may be problematic (see e.g. ^Hl)- In 
addition, the minimal non linear extension of PF theory was found to propagate six degrees of 
freedom rather than five the sixth degree of freedom being ghost like. This increase in the 
number of propagating degrees of freedom was also found to occur in bigravity and in general, 
it may be related to the strong self- interaction of the scalar polarization of the graviton fSl • 

Despite such difficulties, bigravity remains an elegant framework where to investigate the 
massive graviton, particularly because a large class of exact solutions to this theory are known. 
Interestingly, there are solutions where both metrics / and g belong to the Schwarzschild-(Anti) 
de Siter family, which makes them phenomenologically attractive. On the other hand, in dealing 
with a space-time with two metrics, it is natural to ask whether we can make sense of its causal 
structure. In general, the light-cones generated by metrics / and g will not agree, and this may 
lead to pathologies which may restrict the class of physically acceptable solutions. The purpose 
of this paper is to investigate this issue in some detail. 

In Section 2, we review the classification of static spherically symmetric solutions of bigravity. 
These fall into two types. In type I solutions, the area of orbits of spherical symmetry as measured 
by the / metric is 2/3 of the area as measured by the g metric. Otherwise, the solutions are called 
type II. In type II solutions, both metrics can be brought to diagonal form in the same coordinate 
system, whereas this is not always possible in type I solutions, where only one of them can be 
diagonalized at a time. The general solution of type I was given in early work 1^ IT^ . The 
general solution of type II is not known explicitly |22] • We shall consider the subclass where one of 
the metrics is a solution of Einstein's equations with a cosmological constant, for which we discuss 
the general solution. In Section 3, we discuss the causal structure of type I solutions. Section 
4 summarizes our conclusions. Some technical issues concerning type II solutions are left to an 
Appendix. 

^Bigravity (and multigravity in general) is also relevant to the program of "deconstruction" of gravity |11[I12| 
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2 Spherically symmetric solutions in bigravity 

As mentioned in the introduction all type I solutions are known analytically, while the general 
type II solution is not known. In this Section, we review this classification, and comment on a 
subset of type II solutions. 

Following [TU], we consider the action 

—Rn ^ \ f ,A I — z I —R 



2., -.;-y-v-i v^--/ 

Here Lj and Lg denote generic matter Lagrangians coupled to the metrics / and g respectively. 
In solving the field equations, we shall restrict attention to the case where there is only a vacuum 
energy term in each matter sector Lf = —pf,Lg = —pg- There is much freedom in the choice of 
the interaction term in (^. For definiteness, we shall consider the form |10j 



(2) 



with 

u + v = \. (3) 

This reduces to the Pauli-Fierz one in the linear regime, when g is taken to be the Minkowski 
metric (but many other choices could be made with this same property The equations of 
motion, derived from action (^, read 

= ^fTl.. G^,u = (4) 
where G^^p denotes the Einstein tensor, and the energy-momentum tensors are given by 

[^fAr^ - 9''^)(r - gnigaag^r - ga^gar) 

_2(/-/3 - g»''){g^^gf,, - g^^g^p)] + pff^, (5) 

= 2 (f ) [(/"^ ~ g'^^W^ - g''^){ug^ugaagpr - ug^^ga^g^r + ^gf.agaugpr - ^gf^^g^^g^r) 

+2{r^ - g^P) {g^^gp, - g^.g^^)] + pgg^, (6) 



Let us now concentrate on static spherically symmetric solutions. By suitable choice of coor- 
dinates, the most general static and spherically symmetric ansatz for bigravity can be cast into 
the form 

^^dx^dx^ = Cdt^ - 2Ddtdr - Adr^ - B (d^^ + 31^2 q ^^2^ , (8) 

where the metric coefficients are functions of r. As noted in Ref. [21], with the present ansatz we 
have 

DG{, + CGl = 0. (9) 
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Hence, from the equations of motion, the following combination of components of the energy 
momentum tensor must also vanish 



DTi + CTl = ( ^Y{3B - 2r^) = 0. (10) 



B = ^r\ ill] 



Here, we have introduced A = AC + j2Tj. There are two different ways of satisfying this 
equality. Either 

2 
or 

D = 0. (12) 

Solutions obeying the first condition (fTTj) are called type I, and they were first discussed in pUj . 
The remaining ones are called type II and, from (fT^ . they are necessarily diagonal. Let us discuss 
both types of solutions in turn. 

2.1 Type I solutions 

In the case when the vacuum energies of matter fields vanish, the general solution of type I was 
given in pnil21j . The same method can be used in the case when pj and pg are nonvanishing, and 
one readily finds that the general solution is of the form 

g^.Ax^dx" = {l-q) dt^ - (1 - g)-Mr2 - r'^^dO'^ + sin^ ed(p'^) (13) 

/^^dx^dx^ = ^{l- p)dt^ - 2DdMr - Adr^ - \r^{de^ + sin^ ^d(/)^), (14) 

where 

A = ^{l-qr\p + P-q-l3q), (15) 

^'={j^ (i-g)"'(p-g)(p + /9-i-/5g)- (16) 

Here /3 > is an arbitrary constant, and the potentials p and q are functions of r. Substituting 
this ansatz into the expresions (jSHSl); it can be checked that the energy momentum tensors take 
the form of cosmological terms: T^^, = (A//k/)/^j,, and T^^ = {Ag/ Kg)gfj_i^, where 

Kf 4 V2 



Kg 4 



{3v + 913(1 -v)} + pf, (17) 
(^0 (3-^{3u-9(3{l + u)} + Pg. (18) 



It is clear from the equations of motion (0)) that the metrics / and g must belong to the Schwarzschild- 
(A)dS family. Note that the corresponding cosmological constants ()17II18|) are not determined 
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solely by the vacuum energies pf and Pg. They also contain a contribution from the interaction 
term in the Lagrangian. This contribution depends not only on the parameters ( and u (recall 
that V = 1/2 — u), but also on the arbitrary integration constant j3. 

It is somewhat surprising that the cosmological constants depend on an integration constant. 
This situation is reminiscent of unimodular gravity, where the equations of motion are the traceless 
part of Einstein's equations. In this case, the vacuum energy does not determine the trace of the 
Ricci scalar j2S] and Einstein's equations are recovered after using of the Bianchi identities except 
for an arbitrary integration constant that can chosen at will [21] (see also [23). One difference 
here is that we have two cosmological constants Af and A^, and we can only choose the value of 
one of them at will. 

The metric fll4|) can be put in a more familiar form defining a new time coordinate t by 



dt = _|d( + eo <ir}, (19) 

where eo = ±1 is defined by the sign retained for D from equation ()14|) . namely by 



D = - qrWiP-Q)iP + P-^~PQ)- (20) 

With such a coordinate change, the line element ()14|1 now reads 

f^^dx^dx'' = ^{(1 - p)dP - (1 - p)-^dr^ - r\d9^ + sin^ 9d(P^)}. (21) 

As is clear from the previous discussion, the potentials p and q will be given by the familiar 
Schwarzschild-(A)dS forms 

2Mf 2Af , 

P= - + ^r\ 22 

r 9 

q = + -fr , 23 

where Mj and Mg are two additional integration constants with the interpretation of mass pa- 
rameters. 

It is tempting to conclude that this non-linear "theory of massive gravity" is phenomeno- 
logically sound, since the vacuum solutions of General Relativity with a cosmological term are 
recovered, without a trace of the vDVZ discontinuity. In this sense, the mass term does not seem 
to act as an exponential cut-off at a finite range. Rather, it contributes to the effective cosmolog- 
ical constant, which tends to bend space-time on a lengthscale of the order of the inverse mass 
of the graviton (which is of order th? ~ n CY . On the other hand, this contribution from the 
interaction term can be compensated for by a finely-tuned contribution from the vacuum energy 
of matter fields, and then we can have an asymptotically flat solution with exactly the same form 
as for massless gravity. 

It is therefore of some interest to understand the global structure of such solutions, whose 
analysis we defer to the next Section. The perturbative analysis and the investigation of stability 
of these solutions are left for future work [27] . Let us now turn our attention to diagonal solutions. 

^ See also the related discussion of Ref. [21] 
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2.2 Type II solutions 

Type II solutions are defined as those which do not satisfy (jllll . and therefore must be diagonal. 
Although the general type II solution is not known, some progress can be made by further assuming 
that one of the metrics, say is a solution of 

Gl, = Ag^u, (24) 

for some constant A. In the Appendix, we show that the only Type II solutions in this subclass 
are such that 

fli,y = 19i,u, (25) 

where 7 is a constant (to be determined below). 

Substituting (j2SI) in (jSHHl) yields a diagonal energy-momentum tensor for both the g and the / 
metrics, with respective cosmological constants given by 

A, = -6/€gC7"'(-l + 7)(-l - 2m + 27m)/(272) + Kgp,, (26) 
kf = -6/€yC7"'"(-l +l){l~2v + 2^v)/{2^^) + Kfpf. (27) 

Eq. implies that the Einstein tensors for both metrics are identical, and thus 

A, = 7A/. (28) 

This is an algebraic equation which determines 7 in terms of the parameters in the Lagrangian. 

For instance, in the case where the vacuum energies in both matter Lagrangians are set to zero 
Pg = pf = 0, Eq. becomes 

(7 - 1)(-2m +(-!-/€ + 2m(1 - /t))7 + 2Ku-f^) = 0. (29) 

Here k = Hg/nf. In Ref. the case with u = was considered. In that particular case, 

the only solution of (j^^ is 7 = 1, and therefore A/ = A^ = 0. This theory has solutions with 
f = g, where / is an arbitrary solution of the vacuum Einstein equations and the maximally 
symmetric background can only be flat space. However, for m 7^ 0, Eq. (|29p has additional zeroes 
corresponding to the vanishing of the second factor, and therefore we have additional maximally 
symmetric bi-de Sitter or bi-AdS solutions (besides the flat space solution 7 = 1), even in the case 
where the vacuum energies pf and Pg vanish. 

Solutions satisfying ()24|1 and exist for arbitrary interaction between both metrics [not 
necessarily of the form ((21)]. The reason is that by assuming the ansatz where ffj^u = ig^iu, we 
necessarily have 

2 SSintjf , g] 

with constant A^, and similarly for . All that will change from one theory to the other is the 
form of the equation Ag{'~f) = 7A/(7) which determines 7. Maximally symmetric solutions of 
the form ()25j) in a class of bi-gravity theories have been previously considered in |2j. A similar 
argument can be made for multigravity theories, where one can always find solutions where all 
metrics are proportional to each other (see e.g. [2H])- 
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As in the case of type I, both / and g solve Einstein's equations with a cosmological term. In 
this sense, there is no sign of a vDVZ discontinuity either, and one finds the usual Schwarzschild- 
(A)dS metrics, without modification. Unlike Type II solutions, here the two metrics share the same 
light cones and therefore their causal structure doesn't pose any novelties. The phenomenology 
and perturbative analysis of such solutions will be discussed in a separate publication j2Z|. ^ 

3 Global structure of Type I solutions 

In dealing with spacetimes with two different metrics, it is natural to worry about their compati- 
bility from the causal point of view. For instance, there is a Type I solution where / is de Sitter, 
while g is Schwarzschild. The conformal diagrams of / and g look rather different, and the question 
is whether they can be combined at all. The purpose of this section is to investigate such com- 
bined causal diagrams, and the peculiarities which they may introduce. We shall see, for instance, 
that when the solution above is maximally extended, one can send a signal from a geodesically 
complete de Sitter space to another geodesically complete de Sitter space, by following geodesies 
of the companion Schwarzschild metric. Another peculiarity is that the maximal extensions we 
shall consider tend to raise issues of global hyperbolicity, even in cases where both metrics / and 
g are separately globally hyperbolic. Nevertheless, we find that there are no blatant violations of 
causality, in the sense that it is not possible to construct closed future directed time-like curves 
by combining geodesies of both metrics. In what follows, we shall illustrate the construction of 
the diagrams through three different examples. 

The basic procedure is to maximally extend the geodesies of each metric, even when they 
reach the conformal boundary of the companion metric. The resulting causal structure will be 
illustrated by representing light-cones of one of the metrics in the conformal diagram of the other. 
The relevance of such light-cones is two-fold. First of all, matter which is coupled to one of the 
metrics must follow trajectories which are always inside of the future light cone from a given point. 
Second, in the limit of very short wavelength, the interaction term between both metrics becomes 
negligible compared with the graviton kinetic term. Hence, high frequency gravitons will follow 
null geodesies of the corresponding metric. 

Before we proceed with the specific examples, we should briefly comment on a particular type 
of " singularity" which arises in certain bigravity solutions (again, even in cases where both metrics 
are separately smooth). Note that the metric (jl4|) becomes complex in regions where < 0. As 
noted in [21], the coordinate singularity at D = can be removed by a change of variables. This 
is of course true, since / is in the family of Schwarzschild- (A) dS metrics, which are everywhere 
smooth (except perhaps at r = when M/ ^ 0). However, it does not seem to be possible to 

^Assuming that H24|l is valid, these are the only solutions of type II. This does not contradict the perturbative 
approach developed by Aragone and Chela-Flores [221 ' ^^'^ considered spherically symmetric and diagonal solutions 
for / taking g as a fixed flat metric. The difference with the present case is that in their analysis g is not dynamical, 
and therefore there is no need for the tensor Tg to be proportional to the metric [as we have assumed, on account 
of Q and H24() ]. The analysis presented in the Appendix does not exclude the existence of diagonal spherically 
symmetric bigravity solutions, where / would have an asymptotic behaviour along the lines of j22j . It only shows 
that if such solutions exist, then the metric g will not be proportional to /; and neither one will be a solution of 
Einstein's equations with a cosmological constant. 
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find a change of variables wliicli would remove the singularity from both metrics at once, in the 
vicinity of the point at which D"^ changes sign, and which would make both metrics real. The 
reason is that there are geodesies of g which invade the regions D"^ < (with arbitrary slope, in 
fact). On such geodesies, the line element with respect to / is generically complex, and since the 
line element is a scalar, this fact cannot be changed by a coordinate transformation. To avoid a 
complex metric, we could try matching type I solutions with type II solutions at D = 0. This 
possibility is currently under investigation ^27] . 

For most of this Section we shall assume = 1, which ensures positivity of D for all choices 
of the potentials p and g, and therefore seems to be the most natural choice [21 . For certain 
potentials, however, there may be other special values of (3 for which the metric is everywhere real 
(as we shall see in Section 3.3). 



3.1 de Sitter with Minkowski 

Let us choose parameters in (jl7lll8|) so that = and Aj > 0. Then there is a Type I solution 
where g is Minkowski and / is de Sitter. The corresponding potentials in Eqs. (|13lll4j) are given 
by 

p = ^ ^2 ^ H^r\ q = 0. (30) 

Note that each of the spacetimes, characterized respectively by the metrics (fT^ and (fT^ with 
the above defined potentials, has a maximal extension which is geodesically complete (trivial in 
the case of Minkowski). However, combining both together will be non-trivial because the static 
coordinates (t, r) (where we also include implicitly the angular part) cover the whole of Minkowski 
space, but not the whole of de Sitter. Hence, the conformal diagram for the extended de Sitter 
space accommodates all points for which the metric g is defined, but the converse is not true. 
To illustrate the causal structure, let us represent the light-cones of metric g in the conformal 
diagram of /. To this end, it is convenient to use Kruskal-type coordinates, (see e.g. [29j ) 

Note that this involves t (and not t), the temporal coordinate in which / is diagonal. Eq. (p?T|) 
maps the interior of the de Sitter horizon Hr < 1 into the quadrant U < 0, V > of the plane 
(f/, V). The future event horizon for an observer at r = corresponds to f/ = 0, whereas the past 
event horizon corresponds to = (see Fig. 1). The quadrant f/ > 0, V > which lies beyond 
the future event horizon, is similarly covered by the change of coordinates 

The remaining quadrants can be obtained by changing the sign in the right hand side of Eqs. 
(|3HI32j) . As usual, we may perform the conformal re-scaling T = arctanh V + arctanh f/, and 
R = arctanh V — arctanh U, so that the in the new coordinates the four quadrants lie in a square 
of finite size (see Fig. 1). The vertical boundaries correspond to r = 0, while the past and future 
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boundaries of the diagram correspond to r = +00 (which is a spacehke boundary). Note further 
that the coordinate system {t, r) only covers the V > corner of the maximally extended de Sitter 
spacetime but also that it accomodates positive and negative values of U, so that it goes beyond 
the future event horizon. Thus, this coordinate system is similar, as far as the de Sitter metric is 
concerned, to the Eddington-Finkelstein coordinates of a black hole. At this point one might worry 
about a possible singularity due to the presence of the horizon. Indeed, as we discussed above, a 
coordinate singularity in one of the two metric cannot always be removed by a coordinate change 
that renders both metrics non singular. Here the situation is different, and in the coordinates 
(t, r), both metrics are smooth and regular everywhere where t and r take finite values. So the 
U = part of the de Sitter horizon in the V > corner does not result in a singularity in the 
bimetric theory. Things are however more involved for the V = part of the horizon, as we will 
now see. 

To this end, let us consider the light-cones in the Minkowski metric. Radial null geodesies are 
simply given by 

t = er + k (33) 

where e = ±1 corresponds to future and past directed null rays respectively. For e = we obtain 
the space-like t = k slices. In order to represent such geodesies in the conformal diagram for 
metric /, let us first express them in terms of i. For the potentials (jHUI) . Eq. (fTUI) reads 

fir 

di = p-'/^t + — ^(/3 - 1 + H'r^/' (3-'/'dr. (34) 
1 — ii^r"' 

For P = I this yields 

1 , 1-Hr 

t = t — r In 

2H 1 + Hr 



(35) 

The integration constant has been chosen so that t = t at r = 0. For (3^1, Eq. ()34|) can also 
be integrated, but the expressions are a bit more cumbersome and we shall omit them in what 
follows. Note that the change of variables (j35p is discontinuous at the de Sitter horizon. This is 
just as well, since the coordinates t,r become singular at r = rj^ = H~^, and we need to consider 
the Kruskal-type coordinates anyway. Substituting in or in (j^ . we have 

U = ( ^^^] e-^(*-^) , V = e^(*-^) . (36) 



Hr + 1 

As noted above, these expressions are valid both for U < and U > (with V > 0), and so they 
cover both quadrants (jHTj) and (jH^ at once. Now, the radial geodesies are easily given in the U, V 
chart (as a curve parametrized by r) by substituting (fH^ into ^^i. 



.Hr + 1 

Future directed null rays of the Minkowski metric t = r + k, are simply straight lines at 45 degrees 

V = e^'' = const. 
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On the other hand, past directed null geodesies e = — 1, as well as the spacelike geodesies e = 0, 
have a rather non-trivial behavior which is illustrated in Fig. ^ For Hr -C 1, the light-cone 
emanating from r = t = (i.e. = 0) has the same shape as in Minkowski space. However, 
at Hr ~ 1 the past directed light-cone opens up and turns around in the f/, V plane. Beyond 
this turning point, "past directed" null rays of Minkowski start progressing towards the future 
in the de Sitter diagram! In particular, at large affine parameter, Hr — > cxd, both space-like and 
past directed null geodesies of Minkowski meet at the upper left corner of the conformal diagram, 
U — > +CXD, — > 0, which belongs to the future boundary of de Sitter. In fact, the future timelike 
infinity of Minkowski is mapped into the upper right corner of the de Sitter diagram, the 
future null infinity 1^ of Minkowski is mapped into the future null infinity of de Sitter (which 
is spacelike), the spacelike infinity and null past infinity 1~ of Minkowski are both mapped to 
the upper left corner of the de Sitter diagram (see figure |21). The situation is more complicated 
for the past timelike infinity i~ of Minkowski. The latter is split into three pieces: a particle 
moving back in time along a r = constant geodesic of Minkowski space-time would either go to 
the upper left corner of the de Sitter diagram if r > r//, to the lower right corner if r < r//, or 
to the f/ = 0, \^ = central point if r = th- However, a given timelike trajectory in Minkowski, 
stemming from the infinite past {t = — oo, r = rj^) can emanate in the de Sitter diagram from any 
point along the diagonal ^ = 0. The latter diagonal is then representing the whole of the past 
r = th infinity of Minkowski. This can be better seen, plotting the null geodesies of de Sitter into 
a conformal diagram for Minkowski. Inverting ()36p. 



t = r + H-^\TiV, r= ^ ^ , , (38) 

H{i-uvy ^ ' 

outgoing (or incoming) null curves are given parametrically in terms of f/ (or F) by taking V = k 
(or U = k). These are represented in Fig. |21 In particular, one sees that past directed U = 
constant null lines can intersect the V = curve anywhere, while they all asymptote the r = r^j 
curve in the Minkowski diagram as t goes to — oo. 

We may then ask whether it is possible to construct a closed time-like curve by combining 
signals which propagate in the / metric with those propagating in the g metric. We defer this 
discussion to Section 3.4, where we show that this is not possible for general type I solutions. 

A similar analysis can be performed for other values of (3. For (3 > 1, D is everywhere real and 
the causal structure is quite similar to the one described above. A minor difference is that the 
light-cones of Minkowski geodesies are not at 45 degrees near the origin (as they were in Fig. 1). 
This can be easily seen from Eq. (fTI?|) . On the other hand, for (3 < 1 the metric becomes complex 
in the region H^r^ < 1 — (3 (see Fig. 

Let us now consider the issue of global structure. As was stressed above, the coordinates 
(r, t) cover the full Minkowski space corresponding to the metric g, but only half of the conformal 
diagram for the extended de Sitter metric, corresponding to V > (see Fig. la). This portion 
is by itself globally hyperbolic, since the t = k surfaces are Cauchy surfaces for all geodesies of 
both metrics in this region. However, the region > is not geodesically complete, since the 
the null geodesies U = const, of de Sitter reach = at finite affine parameter. To obtain a 
geodesically complete space-time, we can match the solution in the upper half of the conformal 
diagram with a solution in the lower half of the diagram. For this purpose we introduce a second 
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(a) 



Figure 1: Causal diagrams when the / metric is de Sitter (left diagram) while the g metric is Minkowski 
(right diagram) and /3 = 1. The dashed curly vertical line of the left diagram represents a sphere of 
constant radial coordinate r. The solid curly vertical line of the right diagram represents the de Sitter 
horizon r = rn plotted in the Minkowski space-time. We also plotted three radial geodesies of Minkowksi 
space-time emanating from the origin r = at i = 0: the thick dashed (blue) curve is a future-directed 
radial null ray from the origin (notice it is also a null geodesic {V = constant) of the de Sitter space-time), 
the thin solid (green) curve with two arrows is a f = radial geodesic, the thin dashed (red) curve is a 
past-directed null ray from the origin. The last two curves are radial geodesies of Minkowski space-time 
but not of de Sitter space-time. The whole of the Minkowski space-time is mapped onto the half of the 
de Sitter diagram verifying V > Note that the past directed null geodesies of Minkowski turn around 
and start moving towards the future boundary of de Sitter space. This behaviour, however, does not lead 
to closed time-like curves, as discussed in Section ITU 

Minkowski space, with metric g' , which will be covered with coordinates r' and t' . The change 
of variables and (jH^ with the substitutions t —t', U — > —U, V —V, maps the full 
range of the coordinates r',t' into the lower half of the de Sitter conformal diagram, below the 
diagonal ^ = 0. The full diagram, represented in Fig. 3, is now geodesically complete. In doing 
such an extension, we mean we are gluing together one Minkowski spacetime to the other along 
the past infinity of the r = th sphere of the former to the future infinity of the r = th sphere of 
the latter. These infinities do not belong to the Minkowski spacetimes, but to their boundaries, 
while they are located in the interior of the de Sitter spacetime. This provides indeed a perfectly 
fine geometric maximal extension, where all geodesies are complete. 

We should add, however, that a maximal extension is usually required to satisfy the equations 
of motion. The bigravity equations of motion are certainly satisfied everywhere in regions I, II, III 
and IV of Fig. 3, but it is unclear in which sense they are satisfied along the diagonal = 0. The 
problem is precisely that we are joining two Minkowski spacetimes [(a) and (c) of Fig. 3] at a locus 
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Figure 2: Causal diagram for de Sitter with Minkowslti, for (3=1. The left diagram is for de Sitter with 
horizon radius rn, while the right diagram is for Minkowski. The dashed thin lines (with no arrows) are 
t = constant lines. The dashed thick line with one (resp. two) arrow is an r = constant curve, with 
r < th (resp. r > th)- The thin solid line with three arrows represents the trajectory of an observer 
sitting at constant radius r = r// in Minkowski spacetime. The thick solid lines with arrows are past 
directed null geodesies of de Sitter space time U = constant curves. The mapping of the infinities (null, 
spacelike, timelike) of Minkowski spacetimes (i^'*^, X^) has been indicated on the de Sitter diagram. One 
of the stricking feature of those diagrams, is that the past time-like infinity of Minkowski is split between 
the upper left corner (for r > rn), the lower right corner (for r < rn) and the diagonal {r = th) of the 
de Sitter space-time. 

which lies at their conformal boundary. It is conceivable that promoting our maximal extension 
to a solution of the equations of motion might necessitate additional input, such as the inclusion 
of some source at the time-like infinity of Minkowski. Furthermore, as we shall see below, the 
extensions are not unique. This may seem surprising at first sight, but as we shall comment later 
on, a similar ambiguity is present in usual General Relativity when a metric must be continued 
beyond a Cauchy horizon. 

Note that the extended diagram. Fig. |31 is not globally hyperbolic. The t = k surfaces of the 
region > are no longer Cauchy surfaces for the whole space-time, since they do not intersect 
causal geodesies in the lower half of the diagram. A surface which intersects all causal geodesies 
should cut through both regions, V > Q as well as \^ < 0. One such surface is, for instance, the 
horizontal line U = V . The problem is that, as can be seen in Fig. 01 there are geodesies which 
intersect this surface twice (such as the past directed null rays from r = t = 0). A formal proof 
that the maximally extended diagram of Fig. 0] is not globally hyperbolic runs as follows. Let us 
restrict attention to radial geodesies. A Cauchy surface must intersect all causal geodesies once 
and only once. Let us assume that such a surface S exists. In particular, S must intersect the 
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(a) 




Figure 3: Causal diagrams when the / metric is de Sitter (left diagram) while the g metric is Minkowski 
(right diagram) and /? = 1/6. Thick dashed (blue) curve, thin dashed (red) curve, and thin solid (green) 
curve with two arrows, are respectively null (for the two first) and spacelike (for the last) radial geodesies 
of Minkowski space-time. The dashed curly vertical line in both diagram is an r = constant curve which 
is the boundary of the region where one of the metrics becomes complex. 



null geodesic = of de Sitter space. By continuity, it will also intersect the null geodesies 
V = const., in the range —S < V < 6, where S is an arbitrarily small positive number. Let us now 
consider the null geodesic of Minkowski space, parametrized by r in Eq. (jHTf) . and let us choose 
the constant k < H~^lnS. It is clear that the incoming radial geodesic (with e = —1) will start 
at the upper left corner of the de Sitter diagram (at r oo), and work its way down towards the 
right boundary of the diagram (at r=0), while V will always remain in the interval < V < 6). 
Hence, the incoming null geodesic must intersect S at least once before it reaches r = 0. At r = 
it bounces and becomes the outgoing null geodesic V = e^^ < 5, which will intersect S once more 
before it reaches null future infinity. Hence, there are geodesies of Minkowski which intersect S 
twice, which simply means that this is not a good Cauchy surface for all geodesies in the extended 
diagram. 

Let us compare the present situation to that in usual GR. As mentioned above, Cauchy horizons 
are also present in certain maximally extended solutions of GR, such as Reissner-Nordstrom or 
Anti de Sitter space. Whenever there is such a horizon, the equations of motion do not suffice to 
continue the solution past it, and we need additional input. Usually, analytic continuation is used, 
or else some boundary conditions at certain time-like boundaries of spacetime are introduced. As 
mentioned above, in the present context it is not clear whether the equations of motion are satisfied 
or not at the Cauchy horizon of the maximally extended solution, but this is precisely because 
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(b) (c) 

Figure 4: Diagram showing the extension proposed in the text for the de Sitter /Minkowski solution. 
Notations are the same as in Fig. ^ By using a second Minkowski space, we can extend the de Sitter 
diagram of Fig. ^ represented by region I and II above, to the lower half, represented by region III and 
IV above. The de Sitter space-time is now geodesically complete, however the whole space-time it is not 
globally hyperbolic, when both metric are considered on the same footing. If we draw a Cauchy surface 
for all the de Sitter geodesies [such as a horizontal line cutting accross the diagram (6)], this surface will 
intersect some of the Minkowski geodesies twice, while it will fail to intersect some others. 

this horizon corresponds to a point in the conformal boundary of one of the metrics. In this sense, 
the situation is no worse than in GR, where we have to prescribe data on certain boundaries in 
order to determine the maximal extension. Another point to consider is that, physically, Cauchy 
horizons tend to be unstable to perturbations, because of large blueshift effects. The same is 
expected to happen in the present context. Note, e.g., from Fig. |1[ that all future directed null 
geodesies of Minkowski in regions III and IV tend to pile up near the Cauchy horizon at = 0, 
suggesting that there will be a large backreaction near that surface once we include perturbations 

Another interesting fact of the bi-metric solution is that the concepts of causal past and future 
are "broadened" , since signals can be transmitted by matter coupled to both metrics. For instance, 
the observers at r = 0, with ^ > can see signals emitted by all other observers, and hence they 
have no future event horizon. Likewise, observers at r = 0, with V < 0, can emit signals which will 
eventually reach all other observers, and hence they have no past event horizon. It is tempting to 
speculate that cosmological bi-gravity solutions, if they can be made sense of, could in principle 
be relevant to the horizon problem. 



(a) 
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3.2 de Sitter with Schwarzschild 

Let us now replace the Minkowski metric by the Schwarzschild one. In this case, the potentials of 
the Type I solution are given by 

p = H\^, q = . 39 

r 

(40) 

Both metrics have now horizon singularities whenever p = 1 and q = 1, corresponding respectively 
to r = rn and r = rs = 2M. Those are coordinate singularities from the point of view of each 
metric considered separately from the other. However, one might be concerned by the possibility 
to remove such singularities from both metrics at the same time. To study this issue, we first keep 
p and q unspecified, and note that the coordinate change (fT^ reads (with /? = 1, which we shall 
assume in the following) ^ 

dt = dt - dr* + dr , (41) 
r* and r* defining "tortoise" coordinates associated with metric / and g respectively by 

dr 

dr* = (42) 
1-q 

dr 

dr = -— . (43) 
I — p 

Thus, introducing the null coordinates v = t — r*,u = t + r* for the metric g, and v = i — r*,u = 
i + r*, for the metric /, one has from the above expression (|41|) 

dv = dv. (44) 

This means that v is null for both metrics, but also that [v, r, 6, (p) are Eddington-Finkelstein 
coordinates for both metric. In such a coordinates system none of the metric is singular at the 
horizons. 

Coming back to the explicit expressions for p and q (jH^ and substituting those in ()19|) we find 

1 /^^ , ^iH^r^-2M)iH^r^ + if3-l)r-2(3M) ^ \ .... 
^' = 7^^' + (r-2A.f)(l-gV) ' 



For /3 = 1, we have 



' = '-'■•-277'" 



1 - Hr 



1 + Hr 

This matches equation ()4H) where, the Schwarzschild "tortoise" coordinate reads 



(46) 



r + 2Mln|l-r/2M|. (47) 



^we only discuss here the case = +1, the other case, which corresponds to a change in the sign of time, 
follows similarly 
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The analog of Eq. (jHBj) is now 




Hit-r*) 



V = e 



H{t-r') 



(48) 



which, again, is vahd both for f/ > and f/ < (with > 0), covering both quadrants 
and of de Sitter, that is to say the region covered by the Eddington-Finkelstein coordinates 
(f , r, 9, 0). The null and spacelike radial geodesies of Schwarzschild can be written as 



this being obviously valid in the whole region covered by coordinates {v, r, 9, (p). In the U, V chart 
these geodesies are given by 



Again, we find that the null geodesies t = r* correspond to V = const., (or v = const) so \^ is a 
null coordinate both in Schwarzschild and in de Sitter. The other radial geodesies, with e = —1, 
have a more complicated form, which is qualitatively represented in Fig. [7| Note that for this 
figure, we have assumed that the Schwarzschild radius rs is smaller that the de Sitter horizon 
radius rn- 

As we discussed previously, and is manifest from Fig. IHl half of the de Sitter diagram (above the 
diagonal) is mapped onto half of the Schwarzschild diagram (below the diagonal), corresponding 
to the region mapped by the Eddington-Finkelstein coordinates {v, r, 9, 0). Both half-diagrams are 
geodesically incomplete, since some geodesies reach the horizons (which dissects the diagrams in 
two) at finite affine parameter. These geodesies can of course be extended by adding new regions 
of spacetime. If one adds de Sitter and Schwarzschild regions, one obtains a "stair-case" diagram 
with an infinite chain of de Sitter and Schwarzschild space-times, two adjacent de Sitter (resp. 
Schwarzschild) space-times being linked together by a common Schwarzschild (resp. de Sitter) 
space-time. Needless to say, there is also a tension in this case between geodesic completeness and 
global hyperbolicity, as we found in the Minkowski-de Sitter case. 

This seems to apply to more general situations where one of the metrics has an horizon which 
is not shared by the other one ^7\. As noted previously, the new metric (new "step") which can 
be added to the stair does not necessarily correspond to the same solution as the one of the last 
step of the stair, since one of the two metrics does not determine uniquely the form of the other. 
Thus, in general we can construct "stair-case" diagrams with steps having different forms. Note 
further, that in the case considered here, the stairs can always be finished by adding a Minkowski 
spacetime, linked to a Schwarzschild space-time along a sphere of radius vh at time-like infinity. 

3.3 de Sitter with de Sitter 

When both metrics are de Sitter, the potentials are given by 



t = er* + k, 



(49) 





(50) 




q = ny. 



(51) 
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For /3 = 1, the analysis proceeds along the same lines as in the previous Subsection, with the only 
difference that the (de Sitter) tortoise coordinate is now given by 



2Ho 



In 



l-H2r 
l + H2r 



(52) 



The corresponding causal diagram is represented in Fig. El 

Aside from the choice (3 = 1, the de Sitter de Sitter solution allows for another way of having 
> for the entire range of r. Indeed, it is enough to have Hf > jSH"^ and /9 > 1 or iff < /5ii| 

and (3 <1. Choosing for example (3 given by 



(3 



Hi 

Hr 



we have 



Hit = H2t- - In 



or Hi{i — f*) + ln(l + Hir) = H2{t — r*) + ln(l + H2r). Thus, the Kruskal coordinates 
the metric p can be expressed in terms of coordinates t and r as 



U 



Hir - 1 



-^H-iit-r*) 



V 



H2r + 1 



,+H2{t-r*) 



H2r + lJ ' \Hir + l^ 

where r* is given by (j52|) . The corresponding diagram is given in Fig. ^1 



(53) 

(54) 
for 

(55) 



3.4 Closed time-like curves? 

An interesting question regarding the bigravity solutions is whether we can construct closed time- 
like curves by patching together geodesies corresponding to both metrics. 

For /3 = 1 it is easy to show that this cannot be done, by using the "tortoise" coordinates r* 
and f* that we defined in equations (j42j) and (j43|) . as well as the null (for both metric) coordinate 
V (in all this subsection, we keep the functions p and q unspecified). The radial null and time-like 
geodesies of both metrics are given by 

t = er* + k, i = er* + k, 

(Here e = ±1,0 for outgoing and incoming null rays, or for spacelike geodesies, respectively, and 
similarly for e). Thus, any future directed causal curve with respect to f or g has the property 
that df > 0, and dv vanishes only along the outgoing null radial geodesic. Once v increases, even 
if it is by just a little bit, it is impossible to go back to the original value by following a future 
directed time-like curve, which means that such curve cannot be closed. 

Here, we disregard the possibility of making global identifications in the coordinate v, which 
might allow for the construction of a closed loop. Of course, even in flat space with a single 
metric, closed time-like curves could be constructed by global identifications, and in what follows 
we shall ignore this somewhat artificial setup. We shall only be concerned with the possibility of 
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locally constructing closed time-like curves within a given coordinate patch of space-time, without 
identifications. 

To analyse the general case /? 7^ 1 it is convenient to separately consider the following regions 
of space-time: 

a: For [1 — p) < 0, and (1 — g) < the condition dr = defines a space-like surface for both 
metrics / and g. This means that r can only change monotonically along time-like curves of both 
metrics, making it impossible to close them in this region. 

b: For (1 — p) < and (1 — g) > 0, the condition dt = defines a space-like surface for the 
metric g. Also, from ()19|) with dt = 0, we have 



dt 



dr* 



1 + 4 



1 — p 



l-q 



Since t is space-like in metric / this means that the surface dt = [which is also defined by Eq. 
(fH?)|) ] is space-like in metric / too. Hence, t changes monotonically along time-like curves of both 
/ and g, and as a consequence such curves cannot be closed. 

c.- If (1 — > and (1 — g) < 0, then the surface dt = is space- like for /. From (fT^ with 
dt = 0, we have 



dt 



dr* 



1 + /3 



1 



> 1. 



(57) 



1 — p J ' ' \L — p^ 

Since t is space-like in metric g, Eq. (j57|) means that the surface dt = is space-like in metric g 
too, and t must be monotonic on time-like curves, which therefore cannot close. 

d: Finally, if (1 — p) > and (1 — g) > 0, then we must distinguish two cases. For p > g, it 
is easy to see that A > in Eq. (|T^ . and therefore dt = is space-like for both metrics / and 
g. Hence, t is monotonic for time-like curves of both metrics. On the other hand, for p < q, Eq. 
(p77|l for dt = leads to 

dt " 



dr* 



< 1. 



(5^ 



Since now t is time-like in metric g, this means that dt = is a space-like surface for this metric. 
Of course dt = is also space-like for /, and so t is monotonic along causal curves for both metrics. 

This completes the proof for the individual regions listed above. It is remarkable that in spite 
of the strong differences in the light-cone structure of both metrics, it is not possible to draw 
closed time-like curves in any of the regions. The reason is that the future light-cone for one of the 
metrics never contains a part of the past light-cone for the other metric. Thus, we can always find 
a coordinate which labels hypersurfaces which are space-like for both metrics. This coordinate 
must grow monotonically along time-like curves. 

By continuity, at the boundaries in between the regions, the future light-cone of one of the 
metrics can at most touch the past light-cone of the other metric, sharing perhaps a common null 
direction for both metrics. Even if this were the case, a future directed time-like geodesic with 
respect to one of the metrics can never get to the inside of the past light cone with respect to the 
other metric, and closed time-like curves cannot be constructed even if we cross the boundaries 
between the individual regions. 
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In the examples we have examined, the situation where the future hght-cone of one of the metrics marginally 
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4 Conclusions 



In this paper we have analyzed the global structure of a wide class of spherically symmetric 
bigravity solutions. We have reviewed their classification, pointing out that vacuum energy terms 
of the matter lagrangian can be included in type I solutions, and we have discussed a new subclass 
of type II solutions. 

In order to probe the global structure of the solutions, the basic strategy has been to map 
the radial null geodesies of one of the metrics (say g) in the conformal diagram of the companion 
geometry (corresponding to metric /). The relevance of such geodesies is two-fold. First of all, 
they are followed by test particles of the matter sector coupled to each metric. Second, in the 
limit of high frequency, the interaction term between / and g is negligible compared with the 
kinetic term, and short wavelength gravitons will travel along the light cones of the corresponding 
metrics. 

The results illustrate several intriguing features, originating in the different light-cone structure 
of metrics / and g. For instance, a past directed null geodesic of metric g may well end up at the 
future boundary of /, as happens in the case when / is de Sitter and g is Minkowski (Fig. 1). 
In spite of this seemingly bizarre behaviour, we have also shown that for all known solutions, it 
is not possible to construct closed time-like curves by piecing together geodesies of both metrics. 
The basic reason is that the future light cone with respect to one of the metrics never invades the 
interior of the past light cone with respect to the other metric. 

Another peculiarity is that, even if metric / is geodesically complete by itself, some of its 
geodesies reach the conformal boundary of 5^ at a finite value of their affine parameter. In order 
to extend them, further "copies" of the full conformal diagram of g can be incorporated. In 
some cases, this leads to an infinite chain of diagrams, such as the "stair-case" of Fig. |Hl where 
adjacent de Sitter diagrams (corresponding to metric /) are linked together through a common 
Schwarzschild diagram (corresponding to metric g). In this maximally extended bigravity solution, 
one could send a message from one de Sitter space to the next by travelling along a geodesic of 
metric g. 

The inclusion of several copies of the conformal diagram of a given metric leads, in general, to 
lack of global hyperbolicity, even when the "building block" diagrams are globally hyperbolic with 
respect to their own geodesies. Note, however, that global hyperbolicity breaks down in many 
exact solutions of ordinary General Relativity (GR), such as Reissner-Nordstrom or in Anti-de 
Sitter space. This is not necessarily a problem for practical applications, because usually only a 
part of the diagram will be relevant. 

The concept of event horizon is broadened in the context of bi-metric solutions. In principle, 
through gravitational interactions, any observer will be affected at some level by signals propa- 
gating in either metric. In particular, a geodesic observer sitting at r = in de Sitter space will 
have no future event horizon, provided that he or she is sensitive to the signals which are sent 

touches the past hght-cone of the other metric at the boundary between regions does not arise. If it did, then 
there might closed future-directed null curves at such boundary. Note, however, that since the boundary is at 
r = const., this situation can only happen when both metrics have a common event horizon at the same value of 
r. The possibility of having closed null curves on these boundaries may require a case by case analysis, and is left 
for further research. 
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along the geodesies of Minkowski (see Fig. 1). 

In summary, the global structure of bimetric solutions presents some unusual features which, 
nevertheless, do not seem more pathological than those encountered in maximally extended solu- 
tions of GR. It is interesting that all the bi-metric solutions considered are such that both metrics 
/ and g are solutions of Einstein's equations with a cosmological term. Their possible phenomeno- 
logical relevance seems therefore worthy of further consideration, and will be discussed elsewhere 
[?7j (see also ref. [71). 
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Appendix A: Type II solutions 

Here we show that the most general Type II solution where one of the metrics satisfies fj24|) . is 
such that ffj,y = ■yg^u, where 7 is a constant whose value is given by the equations of motion. 

From (j211), we have KT^f- + JT^^ = 0, and plugging expressions (|7j) and (jSJ into Eqs. ©, we 
have 

KTi + JT^^ = y (AJ - CK){W - 2r') = 0. (59) 

Since we are now assuming that B 7^ (2/3)r^, it follows that 

AJ - CK = 0. (60) 
Hence, from (|7|) and (jH)) plugged into ©, 



ATI, + CTl = ^-^j {AJ - CK){3B - 2r') = 0, 



(61) 



and from the equations of motion 

AR{t + CRl = 0. (62) 
From this we obtain (see e.g. pi] for the explicit expressions of the Ricci tensor components), 

/, B'^ A'B' , , 

+ + ^ = °- 

A first integral is given by 

^ = 4a^A (64) 
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where a is the constant of integration. 

Let us now consider the hnear combination 



2Kr 



which again must vanish if is a is a solution of ()24|1 . Thus one either has 

BK + Ar"^ = SAB, (66) 

or 

BJ = Cr^. (67) 

In both cases 

BTl + CT/, = ^ (^4^ J (^^ - Cr^)(i?K - SAB + Ar^) = 0. (68) 

Note that (jHTj) and imply that T^^, = H{r)f^^. The equations of motion require that T'^ 
must be covariantly conserved, which implies that if is a constant. Therefore, / is a solution of 
Einstein's equations with a cosmological constant. 

Consider first the case when ()66|) is satisfied. From this equation and (jUUI) . we can eliminate 
A and C as functions of B and J = . We get from (j64p 



B" Aa' 



(69) 



53 (35-r2)2- 
With the change of variable 

B{r) =r'^F'^{r), (70) 
the differential equation is written as 

rF' = — F (71) 

(3F2-1) ' 

which can be easily integrated to give 



cr = if^^^^)^, (72) 

where c is an integration constant. Notice that 

F(r) = (Vl2 + a2 + a)/6, (73) 
is a solution for a > 0, c — >■ 00 and for a < 0, c = 0, which means 

5 oc (74) 
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In fact, as we shall see, Eq. ([7^) must hold in general. The equation of motion BRf^ + CR^g = 
takes the form [21] 

BC" - CB" + 2 A + {CB' - BC) ^ = 0. (75) 



From ()60p and ()66|) . we have 



and hence 



•4 = #^' C=#^. (76) 

6B — SB — 



B' 



- (35 - r2)2- ^'^'^^ 

Now, Eqs. (j76j) and (j77j) can be used in (j75|) in order to eliminate A and C in terms of B and 
its derivatives (as well as the known function J and its derivatives). The derivatives of B can 
be eliminated from ()64j). and with this Eq. ()75|) becomes an algebraic equation relating B and 
r. Substituting B = r'^F'^, and then eliminating r from Eq. ()72|) . we find an algebraic equation 
involving only F and the integration constants a and c. It turns out that this algebraic equation 
does not vanish identically. Indeed, the first terms in an expansion in powers of F are given by 



(3F(r)2- l)4rVVVl2 + a2 + 

I VVl2 + a2 + a/ VVl2 + a2 + a/ ^ / ^ 

where we have used J = 1 — 2M/r + Agv'^/S. For the zeroth and first order to cancel identically, 
one needs 

A, = 0, M = ^( -glg-° )^. (78) 
-J ^ V 12 + + 

but then going to the next order in F the expression (j75p does not cancel for any value of a. Thus, 
F is fixed to be a constant whose value is determined by ()75|) . From this (f7^ follows. Now, it 
is easy to show that whenever B (x. both metrics must be proportional to each other. Indeed, 
it follows from Eq. (jMjl that A = AC = const, and B = {a'^A)r'^. Also, using JK = 1 and (jUUj) 
we have A = A^^'^K and C = A^^'^J. On the other hand, for constant A, Eq. (|75p reads 

BC" - CB" + 2A = 0. 

Using B = (a2A)r2, C = A^/^J and J = 1 - 2M/r + Agv'^/S, where M and are constants, it 
follows immediately that = A~^^'^, which implies B = A^^'^r'^. It is then clear that f^j^^ = 'yg^^^ 
where 7 = A^/2 jg ^ constant, as we intended to show. 

^^Provided, of course, that the algebraic equation has any solution at all. Otherwise there simply aren't any 
solutions under the assumption H66|l . Note, in particular, from (|72|l and the subsequent discussion, that the 
constancy of F can only be achieved for very special values of the integration constants, but these turn out to be 
the only relevant ones. 
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Next, let us consider the case (jU7j) . Here, we can use (plj) and (jB^ to obtain 



1 

-B^^-v (79) 



2 

B = — TT. (80) 



and equation (fTHjl yields 



Since we have assumed that g satisfies Einstein's equations with a cosmological constant, Eq. 
T^^ should be proportional to g^y with a constant proportionality factor. This is achieved only 
for a = which means C = 7J. This means that both metrics will be proportional, with 

fi^u = IQi^v (81) 

This completes our proof. 

As discussed in the text, the remaining equations of motion determine the constant 7 in terms 
of the parameters in the Lagrangian. 
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(a) (b) 

Figure 5: Same as Fig. 01 witli radial geodesies of de Sitter plotted instead of those of Minkowski. The 
thick dashed (blue) curve is a future-directed radial null ray from the origin (r = 0,t = 0). The thin 
solid (green) curve is a t = radial geodesic of de Sitter. The thin dashed (red) with one arrow curve 
a the past-directed null geodesic from the origin. We also plotted, as thin solid (green) curves with two 
and three arrows, the continuation of the t = curve beyond the horizon r = rn- When mapped into 
the Minkowski diagram, the past directed null geodesies of de Sitter, of region I, reach the timelike past 
infinity of the Minkowski space-time at a finite value of their affine parameter in de Sitter, namely when 
they cross the de Sitter horizon r = rn- Nevertheless, we can "smoothly" continue them in the newly 
added Minkowski solution onto which regions III and IV of de Sitter space-time are mapped. 
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Figure 6: Causal diagrams when the / metric is de Sitter (right) and the g metric is Schwarzschild (left). 
The notations are the same as in figure [21 The main difference with the case depicted in this last figure 
is the presence of the Schwarzschild horizon. The part of the Schwarzschild horizon shown as a thick 
gray line on the right diagram above is mapped to the thick gray line of the left diagram. The part of 
the Schwarzschild horizon which is the diagonal of the right diagram orthogonal to the thick gray line 
is mapped to the upper right corner of the de Sitter diagram in analogy to what was found to happen 
for the de Sitter horizon when the other metric is Minkowski. This shows the possibility to extend the 
Schwarzschild space-time through another de Sitter spacetime joined to the other by the future infinity 
of a r = rs sphere {rs being the Scharzschild horizon) 
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Figure 7: Causal diagrams when the / metric is de Sitter (right) and the g metric is Schwarzschild (left) 
showing the extension proposed in the text for the Schwarzschild space-time. Various radial geodesies 
of Schwarzschild are mapped onto the de Sitter diagram The dashed vertical curly line in the de Sitter 
diagrams indicates the Schwarzschild horizon. Note that we can "send a signal" from region I of the 
lower de Sitter space to region IV of the upper de Sitter space by using the left-moving null geodesic of 
Schwarzschild (thin dashed (red) line). 
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Figure 8: This shows a possible maximal extension of the bi-metric space-times, following the procedure 
given in the text, when one of the metric is de Sitter while the other is Schwarzschild. We are led to 
the "stair-case" diagram, an infinite chain of de Sitter spaces linked to each other through a common 
Schwarzschild diagram. 
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